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The role of elastic interactions between Jahn- Teller ions in formation of various orbital- and charge- 
ordered structures in manganites and related compounds is analyzed. It is shown that such interac- 
tions alone are often sufficient to reproduce the structures observed in different regions of the phase 
diagram. A special attention is focused on stripe structures at high doping levels. 
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I. INTRODUCTION 

As is well established now on many examples, an or- 
bital degeneracy and Jahn- Teller effect play very impor- 
tant role in various transition-metal and rare-earth com- 
pounds In concentrated systems, it leads typically 
to lowering of lattice symmetry, often taking the form 
of a structural phase transition; it is associated with 
(or caused by) the orbital ordering. Correspondingly, 
also magnetic properties are largely determined by the 
type of occupied orbitals according to the Goodenough- 
Kanamori- Anderson rules, see e.g. jl],[|,[| ■ 

Manganites Lai.^M^MnOa (M = Ca, Sr, Ba) with the 
colossal magnetoresistance (CMR) form a very popular- 
nowadays system, for which the orbital effects often play 
very important role. There may be other rare earths (Pr, 
Nd, Sm) instead of La in these series; there exist also lay- 
ered analogues of three-dimensional perovskite mangan- 
ites, often having similar properties. It is established that 
there exists an orbital ordering of some kind in almost all 
parts of the phase diagram of these systems. Thus, it ex- 
ists in Lai_ x Ca :E Mn03 at x = 0; probably in the charge 
ordered insulating phases at small x (x ~ 0.1 — 0.2); at 
the checkerboard-like charge ordered phase at x = 0.5 (so 
called CE phase); in the overdoped regimes 0.5 < x < 1, 
where charge and orbital stripes ||f7| or bistripes || were 
observed. The only exception may be the "optimally 
doped" ferromagnetic metallic phase 0.25 < x < 0.5, 
where the very phenomenon of CMR is taking place: no 
apparent orbital ordering of usual type is observed there 
at low temperatures || although the possibility of an or- 
bital ordering of a new type (with "complex orbitals" ) 
was suggested even for this phase jL0|. In the present 
paper, we will not discuss these more exotic possibilities 
and will limit ourselves to the conventional types of or- 
bital ordering with real orbitals. 

There exist in principle two different mechanisms of or- 
bital ordering, or cooperative Jahn- Teller effect. First, it 
is the usual interaction of degenerate electrons with the 
crystal lattice ^,0,E2[, which is usually considered as a 



source of the Jahn- Teller effect. In concentrated systems, 
this interaction provides the coupling of orbital occupa- 
tion at different sites and can lead to an orbital ordering 
simultaneously with the corresponding lattice distortion. 
The second mechanism is the electronic (exchange) one 
|p^|,p| , which can give rise to both the spin and orbital or- 
derings, and which was quite successful in explaining the 
magnetic and orbital structures in a number of materials, 
notably cuprates |J.Most of the previous theoretical con- 
siderations were confined to undoped systems containing 
integer number of electrons at the transition metal sites: 
d? for Cu 2+ , d A for Mn 3+ . An important class of systems 
is represented by the doped oxides like Lai-^Ca^MnOs, 
where, as mentioned above, similar phenomena of orbital 
ordering are also often observed. One should generalize 
the models used for undoped systems to such cases, in 
particular to be able to explain different superstructures 
observed in these systems, e.g., stripes or paired stripes 
(bistripes). 

Such treatment was recently initiated in where 
we considered the formation of superstructures, includ- 
ing stripes, due to elastic interaction, when we dope our 
system, i.e. substitute ions of one valence (and one size) 
by an "impurity" with the different valence and different 
atomic volume. It was shown there that due to a spe- 
cific nature of these interactions (long-range character, 
~ 1/-R 3 , and, most importantly, different sign - repul- 
sion or attraction - in different directions) the structures 
of different kinds - ID stripes, 2D "sheets" - can be nat- 
urally formed in these cases. 

When considering the case of Jahn- Teller (JT) sys- 
tems, one should modify this treatment by taking into 
account the anisotropic nature of corresponding elec- 
tronic states, i.e. anisotropic electronic charge distri- 
bution(having quadrupolar character) and, correspond- 
ingly, the shape of a Jahn- Teller ion in a respective 
orbital state. Thus, in contrast to the treatment of 
fLlj , we have to consider not the interaction of spherical 
impurities ( "sphere-in-the-hole" model) but rather that 
of anisotropic impurities ( "ellipsoid- in-the- hole"). This 
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generalization was shortly mentioned at the end of [JL4J ; 
in the present paper, we consider this situation in de- 
tail and discuss the applications to the superstructures 
observed in manganites in different regions of the phase 
diagram. As we shall see, in some cases such interactions 
alone are sufficient to reproduce the observed structures. 
In other situations, some extra factors have to be invoked, 
but in any case one can say that the elastic interactions 
considered in this paper are definitely very important in 
stabilizing the observed orbital and, consequently, mag- 
netic structures. 



II. STRAIN MECHANISM OF 
SUPERSTRUCTURE FORMATION 

First, we shortly summarize the results concerning the 
ordering of "spherical" impurities due to elastic forces 
When we put into the host matrix an ion of dif- 
ferent size, e.g. by substituting an ion with the atomic 
volume vq by an impurity with the volume v ^ vq, this 
creates a strain field, which decays rather slowly, as l/R 3 
p5| . Another such impurity interacts with this strain, 
which therefore leads to a certain coupling between im- 
purities. If both the host and impurity ions are spheri- 
cal, this is what is called "sphere-in-the-hole" model. It 
is known that the strain-induced interaction vanishes for 
the isotropic media Jig ] (except of the infinite-range in- 
teraction caused by the "mirror forces" related to the 
sample surface). However, the real crystals are always 
anisotropic. This gives rise to an interaction, which has 
the following form for weakly anisotropic cubic crystals 
0Q: 



V(r, r') 



CQ 1 Q 2 d[ 



R 3 



(1) 



Here Qi — Vi — vq are the "strengths" of the impuri- 
ties, R = |r — r'| is the distance between them, n x , n y , 
and n z are the direction cosines of vector R. The most 
important parameter entering (1) is 



d = Cli — C%2 — 2C44 



(2) 



where Cy are elastic moduli of the crystal. This param- 
eter carries the information about anisotropy: the case 
d = corresponds to an isotropic medium. 

The most significant feature of the interaction (1) is 
that it is attractive in certain directions independent of 
the impurity type (v > vq or v < vq) and of the sign of 
the coefficient d. Thus for d > the interaction is attrac- 
tive along cubic axes [100], [010], and [001] and repulsive 
along face and body diagonals of the cubic cell, [110] and 
[111], and vice versa for d < 0. As argued in uM , this 
quite naturally leads to the formation of superstructures 
(ID stripes, 2D sheets) in insulating systems: the second, 
third, etc. impurities "migrate" toward the first one along 



certain directions, e.g. [001] for d > 0, and finally form 
vertical (diagonal for d < 0) stripes along this direction. 
Important is also that here this motion of impurities is 
provided by the electron hopping and does not require 
real diffusion of atoms; thus, e.g., holes moving in nickc- 
lates transform Ni 2+ ions into Ni 3+ "impurities" , which 
can form stripes due to this mechanism. 

When considering such processes in systems with Jahn- 
Teller ions, e.g. in manganites, one has to generalize this 
treatment to the case of anisotropic impurities - instead 
of "sphere-in-the-hole" we have to consider "ellipsoid-in- 
the-hole" . In this case, there would be two factors deter- 
mining the interaction between such centers: one is the 
elastic anisotropy of the host lattice, and another, spe- 
cific for the JT case, is the dependence of the interaction 
on a relative orientation of the corresponding orbitals. 
In general, one has to consider both factors, but in the 
weakly anisotropic crystals it is often the second factor, 
which plays the dominant role, whereas the first one may 
differ from system to system. 



III. THE MODEL AND RELEVANT 
INTERACTIONS 

We consider below the situation with doubly degen- 
erate e g orbitals, with one electron on them (e.g., ions 
Mn 3 + (4 ej), 



low-spin Ni 3+ (t;L e g)) or with one e g 



As the basis orbitals one can take 
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hole, Cu 2 + {t%e 3 ). 

d z 2 — |3z 2 — r 2 ) and d x -i_ y 2 = \x 2 — y 2 ). The electron 
occupation of these orbitals corresponds to a quadrupolar 
distribution of electron density; the elongated electron 
ellipsoid for the electron z 2 orbitals, and the flattened 
(compressed) ellipsoid for x 2 — y 2 orbitals. One may ex- 
pect also the corresponding displacement of the ligands 
(e. g., the distortion of 06 octahedra): local elongation 
for z 2 electron and x 2 — y 2 hole orbitals, and local com- 
pression in the opposite case. These local distortions will 
lead to a strain field in a crystal, which, as discussed 
above, will provide the mechanism of coupling between 
such ions. The latter can finally lead to the formation of 
one or another type of superstructures. This is the main 
mechanism considered in this paper. 

There are, in principle, several aspects in the electron- 
lattice interactions. The strain-induced interaction men- 
tioned above is the interaction via long-wavelength pho- 
nons. This elastic interaction depends on the type of 
"impurity" and also on the anisotropy of the crystal. For 
anisotropic impurities, like JT centers considered in this 
paper, this interaction depends not only on the relative 
direction between sites in a crystal, but also on the rela- 
tive local distortions on each site. As a result, this inter- 
action becomes very complicated [ ^6|Jl7| ]. In addition to 
long-range interactions mediated by the strain field in a 
crystal, there exists also the interaction mediated by the 
short-range, or optical, vibrations. The detailed form 
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of these interactions was studied in jL2|, and below we 
shall take these interactions into account. It turns out 
that the general structure of both these contributions 
is rather similar: the interactions between the nearest- 
neighbor sites along crystal axes [100], [010], and [001], 
which can be obtained from |L6| and |l2), are exactly 
the same, and the interactions of next nearest neighbors 
along the diagonals [110], [011], etc. have the same sign, 
although somewhat different ratios. Thus, there are two 
contributions, which we consider below, the contribution 
due to the long- wavelength phonons (strain interaction) , 
and that due to the nearest neighbor coupling via short- 
wavelength or optical phonons. 

Generally speaking, one has to include also the elec- 
tronic terms describing the exchange contributions to or- 
bital ordering || and also an electron hopping, which 
can lead to derealization of electrons. We consider be- 
low only the insulating states with localized electrons; 
therefore, we ignore electron hopping here, although in 
principle the tendency to derealization and formation of 
metallic states can compete with the states considered 
below. 

There is one assumption that we make below. In gen- 
eral, an arbitrary combination of the basic orbitals z 2 
and x 2 — y 2 of the type 

\ff) = cos (9/2 |z 2 ) + sin (9/2 |x 2 - y 2 ) 

is allowed, corresponding both to local elongation and 
contraction of Oq octahedra or to any combination of 
these deformations (which would give both prolate and 
oblate ellipsoids or even biaxial ones). Experimentally it 
is well established, however, that only the local elonga- 
tions are realized in practice: out of hundreds known 
Jahn- Teller compounds there are virtually none with 
compressed octahedra (for localized electrons) j|. There 
are also physical reasons for that, related to the lattice 
anharmonicity and higher-order interactions |l^ , [l9tl . In 
accordance to that, we will consider below only the situ- 
ation with locally elongated "impurities" , but the axes 
of local elongations in cubic crystals may be directed 
along x, y, or z axes ([100], [010], or [001]), i.e. angles 
8 may be 0, ±2ir/3 for the one-electron orbitals as in 
Mn 3+ , and 9 = tt, ±7r/3 for the one- hole orbitals at Cu 2+ . 
Thus, we consider the case of elongated impurities (ellip- 
soids) interacting with one another via short-wavelength 
and long- wavelength phonons. The corresponding situa- 
tions are illustrated in Figs, la and lb for the neighbors 
along cubic axes x, y, and z, and in Figs, lc and Id 
- for diagonal neighbors. It is qualitatively clear from 
these figures that the situation of Fig. la (z 2 orbitals at 
nearest-neighbor pair along z axis) would correspond to 
the repulsion (intermediate oxygen 01 feels "conflicting" 
forces), whereas that of Fig. lb would give an attraction. 
This is indeed confirmed by the actual calculations [Q , 
where the constants of these two and of two other pos- 
sible situations were determined, see 1) - 4) in Table 1 



(this is actually the interaction via short- wavelength, or 
optical phonons). 

The situation is somewhat more complicated for diag- 
onal pairs like those in Figs, lc and Id. Here, it is the 
noncentral forces, which determine the resulting interac- 
tions. Still it is qualitatively clear that the interaction 
of oxygens 02 and 03 in Figs, lc and Id would lead to 
an effective repulsion in the first case and to an attrac- 
tion in the second one (note that actually, in contrast to 
our figure, oxygen ions are the biggest ones, so that 02 
and 03 "hard spheres" touch one another, which explain 
our conclusion). The coupling constants corresponding 
to Figs, lc and Id and for two other "diagonal" configu- 
rations (cases 5) - 8) in Table 1) we also take from [ |12[ . 
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TABLE I. The interaction between different orbitals (local 
distortions) for nearest and next nearest neighbors. These in- 
teractions are expressed through two independent constants, 
A and A', characterizing central (cases 1) - 4)) and noncentral 
(cases 5) - 8)) forces. The form and ratios of these interac- 
tions are taken from p"2[ ]; the long-range elastic interactions 
(3) give essentially the same expressions, see text. Here the 
electron orbitals are shown, in which case the local distortions 
have the same form (elongation along the electron cloud) . For 
the one-hole case one should draw the orthogonal hole orbitals 
(e.g. in the case 1) - the hole orbitals x 2 — y 2 ). 



3 



Thus, the resulting nearest-neighbour interactions 1) 
- 8) are expressed through two parameters, A and A', 
corresponding to the central and noncentral forces, re- 
spectively. Note also that these interactions take into ac- 
count both the lattice structure (lattice anisotropy) and 
the anisotropic character of orbital occupation and of the 
on-site distortions. 




(c) (d) 



FIG. 1. Interactions between two Jahn- Teller ions with dif- 
ferent occupied orbitals: (a) and (b) - pair along z axis; (c) 
and (d) - pair along the face diagonal of the perovskite unit 
cell. Displacements of oxygen ions are indicated by arrows. 

The second contribution, due to interaction via the 
long- wavelength phonons for the case of impurities in the 
isotropic medium (i.e., for zero value of parameter d de- 
fined by (2)) can be obtained from the general expression 
presented in pq| ; it has the form 



571-^11 -+- ZC44j-rt ^ 

+2(<®*% +<$><$)} + 

+i^f( 2ff ^-^S-^) (3) 

where en and C44 are elastic moduli, and a aa is a stress 
tensor such that e.g. for the center with the occupied 



orbital 3z 2 — r 2 we have a zz — 1, a xx = a yy = — 1/2 (and 
corresponding expressions for z — > x,y). 

The calculations show that the signs and the ratios 
of the interaction constants for the nearest neighbours 
(cases 1) - 4) in Table 1) obtained from (3), exactly co- 
incide with those via short- wavelength phonons presented 
in Table 1. For the "diagonal" interactions (cases 4) - 8) 
in the Table 1) the signs of the interactions obtained from 
(3) are the same as those given in the Table 1, although 
the ratios of these constants differ somewhat. This is 
probably due to the assumption of the elastic isotropy 
made in (3) (constant d (2) is put = 0), which is appar- 
ently less satisfactory for the noncentral forces. 

For d 7^ 0, the corresponding expression becomes much 
more complicated (see for details). Often, for analy- 
sis of elastic interactions of impurities, further simplifica- 
tions are introduced, namely the longitudinal and trans- 
verse sound velocities are taken to be equal. This leads to 
a rather simple explicit expression for the interaction be- 
tween Jahn- Teller ions with different occupied orbitals, 
situated at arbitrary positions with respect to crystal 
axes |l7]]. For some problems (e.g. for the analysis of 
ESR spectra) , such simplifications work reasonably well. 
However, in our case such approximation is definitely in- 
sufficient, leading to unphysical degeneracy in energies 
of different orbital configurations. At the same time, for- 
mula (3), as it was already mentioned, qualitatively cor- 
rectly reproduces the main features of pair interactions 
even for nearest and diagonal neighbours in the cubic 
crystal. 

Below, in calculating the energies of different ordered 
states, we would take into account only two types of inter- 
actions between the nearest neighbouring JT ions: those 
along cubic axes x, y, z, and along face diagonals. For 
this purpose, we will use the interaction constants pre- 
sented in Table 1 j^J, because they take into account 
both the anisotropy of JT ions and the elastic anisotropy 
of the lattice. 

As we can see from Table 1, the interactions are de- 
termined by two independent constants: the first one, A, 
corresponds to the directions along the cube axes (cases 
1) - 4)) and is determined by central forces between ions, 
and another one, A' - along diagonals of the elementary 
plaquettes - is caused by the noncentral forces [12]. For 
the elastic interactions, these constants are determined 
by different elastic moduli of the crystal: by the bulk 
and shear moduli, respectively. 

In principle, one should include also longer-range in- 
teractions, but their due account would require rather 
extensive numerical calculations, which we postpone for 
future publications. Still, we will show below that even 
using only the nearest neighbor interactions, we can suc- 
cessfully describe many superstructures observed e.g. in 
manganites at different doping levels, and can analyze 
such questions as the relative stability of single vs paired 
stripes, compare the energy of different types of orbital 
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domain walls, etc. One may hope that the longer-range 
interactions, which are still much weaker than those 
taken into account, would not modify our main conclu- 
sions. 



IV. LaMn0 3 , UNDOPED LAYERED 
MANGANITES AND OTHER SIMILAR 
MATERIALS 




First, we consider the simplest case of undoped ox- 
ides, which may be relevant for such systems as undoped 
perovskite LaMnC>3 or KCUF3, layered materials like 
K2CUF4, etc. These systems contain ions Mn 3+ (t^gGg) 
or Cu 2+ (if^e 3 ) at each lattice site, forming a simple cu- 
bic lattice for LaMnC>3 and square perovskite-like lattice 
in layered systems. 

IVa. First, we consider the 2D square lattice relevant 
for layered systems. One immediately sees that if the 



3z 2 



2 ) and 



V 



! ) orbitals are exactly degenerate, 



then with the interactions presented in Table 1 one would 
get the structure shown in Fig. 2. This structure is sta- 
bilized first of all by the attraction (—A/2) of orthogonal 
orbitals (deformations) x 2 and y 2 along x and y direc- 
tions. The diagonal coupling is that of parallel orbitals 
(case 5) in Table 1), i.e. it is also attractive (— A'/3) and 
gives extra stabilization of this configuration. Its energy 
(per site) is 



E = 2(-A/2) + 2(-A'/3) = -A - 2A'/3. 



(4) 



One can easily see that all the other possible orderings 
of these local distortions (locally elongated octahedra) 
give higher energy. The structure described above coin- 
cides with the ordering observed in the layered cuprate 
K2CUF4 [po|l (which was predicted earlier on different 
grounds [|lf). The only difference is that for Cu 2+ such 
a packing of local distortions corresponds not to an al- 
ternation of electron orbitals x 2 and y 2 , but to the corre- 
sponding pattern formed by the hole orbitals (y 2 — z 2 and 
x 2 — z 2 ). The same ordering is also observed in layered 
chromites Ftf^CrCU, see |]3|,p2|. 

Due to a layered structure of corresponding systems 
there may exist in them an initial splitting of the z 2 
and (x 2 — y 2 )-orbitals corresponding to local elongation 
along c-direction, and if this splitting is large enough, it 
can give the ferro-orbital ordering such as the x 2 — y 2 
ordering in La2Cu04. For layered manganites (single- 
layer Lai_ :c Sri +a; Mn04, bilayer L^^Sri^M^Oy), 
we would expect in this case the z 2 occupation at each 
cite. This is apparently the situation in the undoped (x 
= 0) layered manganites p3[. 



FIG. 2. Possible packing of local distortions and orbital 
structure (occupied electron orbitals) for undoped layered sys- 
tems. Oxygen shifts are shown by arrows 

IVb. Let us now consider the case of a cubic per- 
ovskite LaMn03 (we ignore at first the tilting of MnOe 
octahedra leading to the orthorhombicity). One immedi- 
ately sees that the same factors, which acted in layered 
systems, would stabilize the superstructure of Fig. 2 in 
the basal planes of LaMnOs. The complication arises, 
when one considers the ordering pattern in the third (z 
axis) direction. One can have two types of arrangement 
of such planes. They can form either the in-phase struc- 
ture (Fig. 3a), or the out-of-phase structure (Fig. 3b). 
These are usually referred to as the d-type and a-type 
structures, respectively. With the interactions of Table 
1, their energies turn out to be the same: 
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(5) 
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4 



(6) 



Experimentally, both types of ordering are observed 
in perovskites: ordering of the <i-type (Fig. 3a) is found 
e.g. in LaV03, and the a-type structure - in YTi03. In 
KCuF 3 both types of ordering may arise, and stacking 
faults are formed very easily. For LaMn03, only the or- 
dering of the d-type (Fig. 3a) is observed. Apparently, 
there are some other factors, not included on our model, 
which determine the ordering in third direction. 



5 



<<e3 — 




~7 












^ (7 


(a) 





^ 7w\ 

, — — — , 








&~~~\ . 


— e ?' 



(b) 



FIG. 3. Two different kinds of orbital ordering in cubic 
manganites: (a) in-phase and (b) out-of-phase. 

One such factor may be an interplay of the Jahn- Teller 
distortions with the GdFe03-type distortion caused by 
the tilting of the MeC>3 octahedra - one may show 
that for strong enough tilting this interplay stabilizes the 
e?-type ordering of Fig. 3a. In any case, the main motive 
- packing of the distortions shown in Fig. 2 - is the first 
and the most important ingredient in the orbital super- 
structures observed in many perovskites with JT ions, 
and this directly follows from our mechanism. 

Thus, we see that the electron-lattice coupling quite 
naturally gives the correct lattice and orbital superstruc- 
ture of undoped manganites (see also 12 23]) 



V. LOW-DOPED MANGANITES 

There are several interesting systems, which show 
rather puzzling properties at low doping level. Thus, 
Lai-^Sr^MnOs has an insulating ferromagnetic state at 
about x = l/8(~ 0.1 < x < 0.14) Such a state 

appearing in oxides is in general a rare occasion: most of 
the insulating systems are antiferromagnetic, and ferro- 
magnetism usually coexists with (and is explained by) the 
metallicity (double exchange) . Probably, the only way to 
stabilize a ferromagnetic insulating state is by forming a 
special orbital ordering [Q, so as to give ferromagnetic 
exchange according to the aforementioned Goodenough- 
Kanamori- Anderson rules. The superstructure was in- 
deed observed in Lai-^Sr^MnOs in this region pq ], and 
it was interpreted in |^| as predominantly the charge 
ordering - ordering of holes (Mn 4+ ions) of more or less 
bcc type in the Mn 3+ matrix. The remaining Mn 3+ ions, 
of course, should have certain orbital occupation in an 
insulator, which could in principle finally lead to ferro- 
magnetism. 

This problem was considered theoretically in [p9| , 
where it was shown that if one assumes the charge order- 
ing pattern of Yamada et al. p8| ], one would indeed obtain 
certain orbital ordering, which can be treated as an order- 
ing of orbital polarons 1 30 3^] leading to a ferromagnetic 
spin ordering. These orbital polarons are the objects, 
in which the central small-size Mn 4+ ion is surrounded 



by larger Mn 3+ ions, with the lobes of their orbitals di- 
rected toward the central Mn 4+ . The elastic forces are 
definitely rather important in forming such objects and 
in their possible ordering (see also [j52|). However, the 
Hartree-Fock calculations carried out in |2{J] have shown 
that the better structure with the lower energy is the one, 
in which the xy planes alternate as follows: one plane is 
a pure Mn 3+ plane with the superstructure of Fig. 2; the 
next plane contains all the holes (Mn 4+ ions) forming 
"vertical" stripes (stripes along x or y directions), see 
Fig. 4. Such a state turned out to be also ferromagnetic; 
this structure is consistent with the lattice symmetry ob- 
served in p8j . 

Although the numerical calculations of 0] were car- 
ried out using different starting point, one can argue 
that it is the special stability of the orbital ordering 
of Fig. 2 that helps to stabilize (or is even a driving 
force of) such a superstructure in low-doped mangan- 
ites, with segregation of pure undoped Mn 3+ planes, all 
the holes being in other planes. If it is true, one can 
expect similar phenomena also in other systems, e.g. in 
low-doped Pri_ x Ca x Mn03, where the ferromagnetic in- 
sulating phase is also known to exist at 0.1 < x < 0.3 
p3| , [39|] . It is possible in principle that the structure of 
this system, e.g. close to x = 1/4, is also formed by the 
ordered array of orbital polarons |3l[]. But the better 
alternative may be again the segregation of holes in ev- 
ery second plane, undoped planes being indeed like those 
in Fig. 2, and doped planes (which in this case would 
have the effective doping x e s = rihoios = 0.5) having ei- 
ther stripe-like structure or even the superstructure of 
the CE type, characteristic of x = 0.5 (see below). 



Mn 3+ plane 




< doped layer 

(holes in chains) 



Mn 3+ plane 



FIG. 4. A possible orbital ordering at n — 1/8. 

The detailed type of charge and orbital ordering in 
Pri-^Ca^MnOs in the vicinity of x = 1/4 is not known 
yet, although the first measurements by the anomalous 
X-ray scattering |55| have shown that in this system there 
indeed exists a superstructure, which can be consistent 
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with that in Fig. 4 and/or with the CE type superstruc- 
ture of Fig. 5a. It would be very interesting to check this 
possibility experimentally in more detail. 



VI. HALF-DOPED SYSTEMS: CHARGE AND 
ORBITAL ORDERING 

Via. Ground state 

In most of half-doped manganites i^i^^Me^MnOs, 
x = 0.5 (R = La, Pr, . . . ), there exist at low temper- 
atures the charge and orbital ordering accompanie d by 
the antiferromagnetic ordering of the CE type |pi| , p9| , 
the latter usually appearing at still lower temperatures. 
Nowadays, one often uses the term "CE ordering" also to 
denote the corresponding charge and orbital structures. 
This ordering is illustrated in Fig. 5a. It consists of the 
checkerboard arrangement of Mn 3+ and Mn 4+ ions with 
the corresponding orbital ordering at Mn 3+ sites. Of 
course, one should not take notations "3+" and "4+" 
too literally: actual degree of charge disproportionaliza- 
tion can be much less, e.g., ~ (3.5 ±0.2) |37|) . This is the 
structure in the basal plane; the stacking of these planes 
in the third direction is in phase, so that 3+ and 4+ sites, 
as well as the corresponding orbitals, are on top of one 
another. 
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FIG. 5. Charge and orbital orderings in half-doped man- 
ganites: (a) actual ordering (the so called CE structure); (b) 
alternative orbital ordering with all occupied orbitals at Mn 3+ 
of the same type (x 2 orbitals). Empty circles denote Mn 4+ 
ions, the electron density distribution for the corresponding 
orbitals is shown for the Mn 3+ sites. 

There are two questions related to this superstructure. 
The first one concerns the type of charge ordering. It 
is quite natural to expect the checkerboard arrangement 
of Mn 3+ and Mn 4+ ions, or of extra electrons: such or- 
dering is actually analogous to a Wigner crystal, i.e. it 
is favored by the Coulomb forces. These considerations, 
however, are definitely not enough, because in this model 
one would expect charge (valence) alternation also in the 
third direction, which is not the case. The solution of this 
problem may again be connected with the elastic inter- 
actions: it can be shown 0,[l(| that they typically favor 



the inclusion of one phase (e.g. Mn 4+ ) in another (Mn 3+ ) 
in the form of infinitely thin slabs - sheets, oriented in a 
crystal in a particular way, so as to minimize total strain 
energy. These may be just Mn 4+ sheets observed in the 
CE structure. 

Another question concerns the type of orbitals occu- 
pied at Mn 3+ sites. If we assume the checkerboard charge 
ordering, we may have different types of occupied or- 
bitals, e.g., one can in principle get the ordering shown 
in Fig. 5b with all orbitals of the same type (here, x 2 or- 
bitals) instead of alternation of x 2 and y 2 diagonal rows 
( "stripes" ) of Fig. 5a. 

The type of orbital ordering in the checkerboard charge 
structure can be again analyzed using the interactions 
listed in Table 1. Here we also include the interaction 
of two types of nearest neighbors: the first neighbors in 
the Mn 3+ sublattice are here the diagonal interactions 
5) - 8) of Table 1, and the second neighbors are those 
along x and y directions in Fig. 5. The only difference 
is that these interactions are now at twice the distance 
of the nearest neighbors Mn 3+ and Mn + in these direc- 
tions, so that according to the general nature of elastic 
interactions, which decay as 1/-R 3 , these couplings with 
parameter A (cases 1) - 4) in Table 1) should be now 
multiplied by factor 1/2 3 = 1/8. 

From Table 1 one sees that the diagonal interactions 
are attractive for the same orbitals, e.g., x 2 and x 2 (case 
5) in Table 1), but repulsive for x 2 and y 2 (case 6)). From 
this point view, one should rather expect that the orbital 
structure in this case would correspond to Fig. 5b, and 
not to Fig. 5a. However, the interactions along x and y 
directions tend to stabilize the structure of Fig. 5a. One 
can easily calculate the energies of these two competing 
states (per site): 
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(7) 



Here, the first two terms are the diagonal interactions 
with the first neighbouring Mn 3+ ions and the last term 
is the interaction along x and y directions (at a distance 
of two lattice parameters). 

Similarly, the energy of the state presented in Fig. 5b 

is 



E$b — 2 



-A' 
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—A 
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2A' 5A 
--+32 (8) 



Thus, if the central force interaction A is strong enough, 

A'/A < 3/8 = 0.375, (9) 

the CE structure of Fig. 5a would be stable; for the op- 
posite inequality one would get the structure of Fig. 5b. 
Since actually the CE structure is observed in experi- 
ment, one should conclude that inequality (9) is fulfilled 
in manganites (within the approximations made). One 
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can also show that all other feasible types of orbital or- 
derings have higher energy. 

One has to stress that the assumption of the checker- 
board charge ordering is crucial for this conclusion; if we 
would not make this assumption and use only the JT- 
induced interactions of Table 1, one should rather get 
in the model under study the phase separation: all elec- 
trons (Mn 3+ ions) would form a dense cluster with the 
superstructure of Fig. 2. This would minimize the inter- 
action energy of the JT-distorted sites. Thus one has to 
add some extra factors to this model, e.g., the Coulomb 
interaction preventing such phase separation. However, 
assuming the checkerboard charge ordering, we can ex- 
plain the concomitant orbital ordering of Fig. 5a in our 
model. (Nevertheless, the possibility of phase separation 
cannot be completely discarded even in this case - and 
all the more so at x ^ 0.5 0].) 
VIb. Orbital domains 

It is also of interest, and of practical importance, to 
study in our model the "cost" of creation of defects, no- 
tably orbital domain walls. As argued in J34| and espe- 
cially in [ 06| , such domain walls are easily formed in CE 
states, e.g. in Lao.sCao.sMnOa and in Pri_,cCa x Mn03. 

One can easily calculate in our model the energy of 
such defects. The simplest one is the local defect - rota- 
tion of the orbital on one site, e.g. from x 2 to y 2 . With 
our interaction constants it would cost us AE = 9/16A 
(diagonal interactions remain the same for the real CE 
structure corresponding to Fig. 5a). More interesting is 
the situation with the domain walls. There arc four types 
of them, classified and drawn in Figs. 3b - 3e of p6[| . Do- 
mains separated by domain walls 3b and 3c of ]36[ are 
those with the same propagation vector, i.e. the same 
orientation of orbital stripes, and 3d and 3e are domain 
walls between two domains with the propagation vectors 
rotated by 90 degrees (orbital twins). Straightforward 
calculations show that the energies of domain walls (per 
unit length, i.e. per Mn 4+ -site at the center of domain 
wall) are 



AE; 



36 



V + 32 A 



(10b) 



for this coincides with the critical value (A'/A) = 3/8 (9) 
obtained above, which confirms the consistency of our 
treatment. 

An important conclusion is that, at least in our model, 
in the range of existence of the CE structure, (A'/A) < 
(A'/A) c = 3/8, the domain walls of the type 3d and 3e 
are easiest to form, and they have the same creation en- 
ergy. Thus, it will be presumably orbital twins 3d and 
3e which will appear in real systems. The driving force 
for their formation is apparently the same as in ordinary 
martensitic transitions pC| ], i.e. they are elastic domains 
caused by strain in crystals (this is suggested also by A. 
Millis, as cited in |3(|). different from what was usually 
assumed: instead of the 



VII. STRIPES AND BISTRIPES IN OVERDOPED 
MANGANITES 

Finally, we turn to the most controversial problem 
- that of the type and the origin of superstructures 
in overdoped manganites, e.g. stripes and bistripes in 
Lai-zCa^MnOa at x = 2/3; 3/4 First of all, there 

exists an experimental controversy: different experiments 
give for the same system and even for the same samples 
different results. Thus, for x = 2/3 the high-resolution 
electron microscopy gave the results taken as a signa- 
ture of paired stripes or bistripes Q , whereas the neu- 
tron scattering results were interpreted in terms of single 
stripes (called "Wigner crystal" in M). 



Q Q O i O Q ©■ 



O OO O 



|0 $- 

/ 



11 



o — o 



Hf O O OO O 

-e — o oa o — © 

(a) 



Q Q QO Q 



O OO G 



B G-OO 





-e- 





o — $ — 



G 



4 



e — g-oo 



«-ee-o 



© — o Oo G 
(b) 



AE- 



3c 



3 A'4- 9 A 
4 A + 32 A 



AE* d = AE,, = 



A 



M, 



4 32 



(10c) 



(lOd-e) 



As we see from (10), the domain wall 3c is the most 
"expensive" one, it always has positive creation energy. 
The situation with other three types of domain walls, 
3b, 3d, and 3e in the notation of pq |, is different: they 
cost us some energy for small (A'/A), but they can be 
created spontaneously (AE < 0) if diagonal coupling A' is 
sufficiently strong. As we see from (10), the critical value 
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FIG. 6. Possible stripe structures in the basal plane of 
Lai-zCazMnOs at x = 2/3: (a) single stripes; (b) bistripes 
with different orbitals; (c) bistripes with similar orbitals. 
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These two situations, with the corresponding orbital 
ordering, are illustrated in Fig. 6a, b for the case x = 2/3. 
The orbitals at each stripe are parallel, e.g. one stripe is 
formed by x 2 orbitals and another stripe by y 2 orbitals. 
In the paired stripe picture (Fig. 6b), these two diagonal 
stripes, with x 2 and y 2 orbitals, come close together being 
separated by one Mn 4+ diagonal row. 

Again, we have here several questions: why stripes here 
are better than any other type of charge ordering? (e.g., 
for x — 3/4, one can arrange the real "Wigner crystal" 
structure with the Mn 3+ sites forming a face-centered 
2D lattice); and why do stripes have a particular orbital 
structure? And the third, most burning question is which 
state, single or paired stripes, is better? 

One can try to approach these problems using our 
model with the interactions of Table 1, adding to it a 
restriction forbidding the formation of dense Mn 3+ clus- 
ters. Similar to Section VI, we calculate now the energies 
of the single stripes of Fig. 6a, paired stripes of Fig. 6b, 
and also the possible competing type of bistripes shown 
in Fig. 6c, with all orbitals of the same type, e.g. x 2 . 
Using the interaction parameters of Table 1 and keeping 
again the interaction of the first (diagonal) and of the 
second "direct" neighbours at a distance of two lattice 
constants in x and y directions, we obtain: 
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Expression (11) is quite clear. In (12), the first term 
describes the diagonal interaction along the stripe, the 
second - diagonal interactions between stripes in the pair 
(each such bond belongs to two sites, therefore we have 
factor 1/2 in it), and the third term is the next-nearest- 
neighbour interaction of x 2 orbital of one stripe with the 
y 2 orbitals of its "mate" . Similarly, one can understand 
the meaning of different terms in (13). 

Comparing the energies of these three states, one ob- 
tains that there exist two critical values of parameter 
v = A'/A 



,(<**>) - 



= 3/10 = 0.3, 



Xac) _ 



15/32 ~ 0.47. 



(14) 



(15) 



For A'/A < v£ b) = 
are stable. For 0.3 



0.3, the paired stripes of Fig. 6b 



(ab) 



< A'/A < v, 



(ac) 



0.47, the 



single-stripe phase of Fig. 6a would have the lowest en- 
ergy among the three states, which we compare. And fi- 
nally, for A'/A < z^r = 0.47, the paired stripes of Fig. 6c 



would have the lowest energy. We again discard here the 
possibility of phase separation, where all electrons are as- 
sembled in the bulk CE-like phase of Fig. 5a, which can 
arise at small A'/A, and into the phase of the type shown 
in Fig. 5b for large A'/A, assuming that such phase sep- 
aration into electron-enriched and electron-depleted re- 
gions would be prevented, e.g. by the long-range Coulomb 
interaction. 

Comparing these results with those of Section Via 
(half-doped systems), we can see that for realistic val- 
ues of A'/A < 0.375 (9), which are needed to stabilize 
the experimentally observed CE phase, we may have ei- 
ther the paired stripes § (Fig. 6b), or single stripe phase 
of Fig. 6a, see phase diagram in Fig. 7. Unfortunately, 
within this treatment we cannot make a unique choice 
between these two possibilities. Nevertheless, we can say 
that at least the sometimes used arguments ||, which 
ascribe the pairing of stripes to strain interactions, are 
not very convincing: we see that the same strain interac- 
tions for different ranges of parameters may rather favour 
single stripes. Apparently, a more detailed treatment, in- 
cluding also long-range interactions, is needed to answer 
theoretically the question, which state, single or paired 
stripes, would be preferred for real systems (although it 
is, of course, primarily an experimental problem). 



x = 0.5 



x =2/3; 3/4 



FIG. 7. Regions of stability of different phases for x — 0.5 
(upper half-plane) and for overdoped manganites [x = 2/3; 
3/4, lower half-plane). Note that we do not consider fully 
phase-separated states. 

As an indirect argument we may only use the results 
of the study of domain walls at the end of Section VI. As 
shown there, the energy cost of domain wall formation 
goes to zero at the critical point (A'/A) = 0.375 (9). As 
follows from the experimental data ]36[ ], orbital domains 
are formed very easily in half-doped manganites. It may 
be a signature that real systems are not far from this 
critical point. If true, this would mean that, according to 
the results (12) and (13) and to Fig. 7, single stripes may 
be preferable, although the energy of paired stripes may 
be not very far from it, so that they may become stable 
e.g. at the surface of the crystal. This could resolve 
the experimental controversy between the results of |j] 
and 0. 
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VIII. MODELING OF ORBITAL ORDERING BY 
INTERACTING QUADRUPOLES 

As is well known || , the orbital ordering is closely re- 
lated to the ordering of electric quadrupolcs: the asym- 
metric distribution of charge density at a transition metal 
ion with a particular occupation of one of degenerate or- 
bitals corresponds to a quadrupolar charge distribution, 
and, moreover, the interaction of orbitals via lattice dis- 
tortions is very similar to the conventional quadrupole- 
quadrupole (QQ) interaction. It is instructive to see 
whether one can model the elastic interactions listed 
in Table 1 by this QQ interaction. Of course, in real 
transition metal compounds this direct QQ interaction is 
much weaker than that via phonons (it may be relevant 
though for rare earths compounds, where one indeed of- 
ten uses the terminology "quadrupolar ordering" instead 
of our "orbital ordering" and "cooperative Jahn- Teller 
effect"). The distance dependence is also different: the 
QQ interactions decay as 1/R 5 , whereas the interactions 
via phonons - only as 1/R 3 . Nevertheless, the operator 
structure of the corresponding terms in the Hamiltonian 
is very similar, and one may expect that at least the 
signs, and probably the ratios of corresponding coupling 
constants would be similar for both these mechanisms. 

The general form of interaction between two quadru- 
poles Q a p and Q 7< s, located in positions with radius vec- 
tor R connecting them, can be obtained from the general 
treatment [|42fl: 



u QQ ,= 



12R 5 



a/3 



dV 



(16) 



where n a etc. are directional cosines of radius-vector R. 

Using this general expression for axially symmetric 
quadrupoles (Q zz = Q, Q xx = Q yy = —Q/2), one can 
easily calculate the quadrupolar interaction constants for 
the cases 1) - 8) of Table 1. They turn out to be 
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where A = 3Q 2 /2R 5 . In calculating the values of Vn — Vg\ 
in (17), we separated out the factor l/(\/2) 5 , since the 
distance for diagonal pairs is larger by a factor of y/2 
than that for the "direct" nearest neighbours. 

One can see that indeed the signs of the QQ interac- 
tions in (17) coincide with those of elastic interactions 



(see Table 1); this gives an extra support to these val- 
ues used in the main part of our paper. Moreover, even 
the ratios of the interaction constants are rather close, 
at least for cases 1) - 4). In principle, one can try to 
find out the type of the ordered ground state for the QQ 
case, using general expression (16) and taking into ac- 
count the longer-range interactions and interactions for 
all possible directions. This task would require rather ex- 
tensive numerical calculations (or some clever summation 
of the Madelung type); we do not do this here. One can 
only formulate some interesting problems, which should 
be addressed in this field. 

(1) Would there be indeed a long-range ordered state in 
this system? We remind that for a similar problem with 
the dipole, not quadrupole, interactions it was proven 42] 
using the approach similar to that of Mermin- Wagner 1 44 ] 
that there is no long-range ordering for the dipole-dipole 
system in the case of 3D cubic lattice. Such possibility 
is not excluded also for the QQ interaction (for systems 
with quadrupoles at each lattice site). 

(2) Another curious problem is what would be the sit- 
uation with the mobile quadrupoles if their concentra- 
tion is small (their number is less than the number of 
lattice sites). Such quadrupoles definitely attract each 
other at certain directions and for certain relative orien- 
tations of their main axes. Would they form structures 
like stripes, sheets, or some more complicated superstruc- 
tures, and which in particular? Note that this possibility 
is rather similar to the one suggested in jl4| for elas- 
tic interactions of spherical "impurities" in anisotropic 
media. Here, however, one considers the interaction of 
anisotropic centers in the isotropic space. Nevertheless, 
the outcome may be similar (actually the consideration 
of Section IV have much in common with this problem; 
but for the QQ case one can go beyond approximation of 
the nearest-neighbor interactions made above and con- 
sider the problem in its full complexity). Note also that 
similar question exists also for mobile dipoles. 

These situations have not only an academic interest: 
one can imagine that similar superstructures of various 
kinds can arise due to these interactions, e.g. in the so- 
lutions of anisotropic molecules in certain solvents, etc. 

The questions discussed in this Section have only in- 
direct relation to the main problem considered in the 
present paper, but still we see much in common between 
these situations, at least in principle. 



IX. CONCLUSIONS 

In conclusion, we considered in this paper the possi- 
bility of formation of different superstructures in insu- 
lating systems with Jahn- Teller ions like some cuprates, 
chromites, and especially manganites at different doping 
levels. We saw that even keeping only the interactions 
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between a few nearest neighbors, we can successfully de- 
scribe the formation of different orbital superstructures. 
In particular, using only the elastic interactions one can 
immediately obtain the observed ordering pattern of un- 
doped manganites. By assuming the checkerboard charge 
ordering, we can easily explain by this mechanism the 
type of orbital ordering (CE structure) observed in half- 
doped manganites [x = 0.5). Using the same physical ar- 
guments, we suggested that the superstructures observed 
in the low-doped manganites (ferromagnetic insulating 
phases at x = 0.1 — 0.3) involve the segregation of doped 
holes in consecutive planes, this structure being favored 
by the special stability of the orbital ordering of the "un- 
doped LaMn03-type" - the ordering shown in Fig. 2. For 
the overdoped manganites like Lai-^Ca^MnOs, x > 0.5, 
we argued that the stability of stripes or bistripes ob- 
served in this doping range (especially at x — 2/3; 3/4) 
is also determined by the elastic forces, although we were 
not able to determine uniquely which of these alterna- 
tives, single or paired stripes, is more favorable. In our 
approach, the result depends on the ratio of two con- 
stants related to the noncentral and central forces respec- 
tively, which we do not know a priori. But in any case, 
the elastic forces largely determine the very tendency to 
the stripe formation in insulating overdoped manganites, 
as well as their orbital structure. 

As mentioned in the Introduction, there exist also 
other possible mechanisms of orbital ordering, in particu- 
lar, the superexchange interaction PJl3]]. Typically both 
these mechanisms, the exchange interaction and the in- 
teraction via lattice, lead to the same orbital structures 
in cases when we have Jahn- Teller ions at each site. It 
is quite difficult to find the cases, in which the outcomes 
of these two models would be qualitatively different [E5| . 
In some sense, it is gratifying that very simple consider- 
ations presented above immediately give the correct or- 
bital structures for the cases like LaMnC>3 or K2CUF4, 
whereas rather hard work is required to get these struc- 
tures with the electronic (exchange) mechanism. Thus, 
apparently both these mechanisms are operational in real 
systems, and the relative importance of each is still an 
open problem, at least for the cases of dense systems (see 
also jl6| in this context). 

But when we turn to the dilute systems, like doped 
manganites, especially in the overdoped regime x > 0.5, 
then apparently the elastic interactions play the most 
important role, first of all providing the essential mecha- 
nism of stripe formation, and also determining their or- 
bital structure: one can expect that the exchange inter- 
actions, being essentially short-range, are less efficient in 
these cases. It is mainly in these systems, where we be- 
lieve that our approach may be most fruitful. 
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